Observation of the rare 77 — > e + e e + e decay 
with the KLOE experiment 



The KLOE Collaboration 

F. Ambrosino d ' e , A. Antonelli a , M. Antonelli a , F. Archilli 1 ^, 
(Nj '. I. Balwierz b , G. Bencivenni a , C. Bini g ' h , C. Bloise a , 

q ■ S. Bocchetta k ' £ , F. Bossi a , P. Branching, G. Capon a , 

T. Capussela a , F. Ceradini k ^, P. Ciambrone a , E. Czerwinski a 
E. De Lucia a , A. De Santis g ' h , P. De Simone a , G. De Zorzis< h 
A. Denig c , A. Di Domenico §' h , C. Di Donate 6 , B. Di Micco w 



(N 



Oh- 

CD 



X 



5 

M. Dreucci", G. Felici", S. Fiore s ' h , P. Franzini«'\ C. Gatti", 
P. Gauzzi S. Giovannella E. Graziani M. Jacewicz 
J. Lee-Franzini a ' m , M. Martemianov °, M. Martini a ' f ' p , 
P. Massarotti d ' e , S. Meola d ' c , S. Miscetti a , G. Morello a , 
M. Moulson a , S. Miiller c , M. Napolitano d ' e , F. Nguyen k ' £ , 

k,£ 
a 



> . 

■ M. Palutan a , A. Passed , V. Patera a ' f , I. Prado Longhi 

O ! P. Santangelo a , B. Sciascia a , M. Silarski , T. Spadaro 

C. Taccini k '^, L. Tortora^, G. Venanzoni a , R. Versaci a ' f ' q '*, 
o: G. Xu a ' n , J. Zdebik b 



and, as members of the KLOE-2 collaboration: 

D. Babusci a , D. Badoni^, V. Bocci h , A. Budano k ' f , 
S. A. Bulychjev , P. Campana a , E. Dane a , G. De Robertis s , 
D. Domenici a , O. Erriquez r ' s , G. Fanizzi r ' s , F. Gonnella 1 ^, 
F. Happacher a , B. H6istad v , E. Iarocci f ' a , T. Johansson v , 
V. Kulikov , A. Kupsc v , F. Loddo 8 , M. Matsyuk , R. Messi iJ , 
D. Moricciani-i, P. Moskal b , A. Ranieri s , I. Sarra a , 
M. Schioppa^, A. Sciubba f ' a , W. Wislicki w , M. Wolke v 

a Laboratory Nazionali di Frascati delVINFN, Frascati, Italy. 
b Institute of Physics, Jagiellonian University, Krakow, Poland. 
c Institut fur Kernphysik, Johannes Gutenberg - Universitdt Mainz, Germany. 
d Dipartimento di Scienze Fisiche delVUniversitd "Federico IF', Napoli, Italy 



Preprint submitted to Elsevier 



26 January 2013 



°INFN Sezione di Napoli, Napoli, Italy 

f Dipartimento di Scienze di Base ed Applicate per I'Ingegneria dell'Universita 

"Sapienza" , Roma, Italy. 

g Dipartimento di Fisica dell'Universita "Sapienza" , Roma, Italy. 

h INFN Sezione di Roma, Roma, Italy. 

1 Dipartimento di Fisica dell'Universita "Tor Vergata", Roma, Italy. 

3 INFN Sezione di Roma Tor Vergata, Roma, Italy. 

k Dipartimento di Fisica dell'Universita "Roma Tre", Roma, Italy. 

e INFN Sezione di Roma Tre, Roma, Italy. 

m Physics Department, State University of New York at Stony Brook, USA. 

11 Institute of High Energy Physics of Academica Sinica, Beijing, China. 

Institute for Theoretical and Experimental Physics, Moscow, Russia. 

p 'Present Address: Dipartimento di Scienze e Tecnologie Applicate, Universitd 
Guglielmo Marconi, Roma, Italy. 

^Present Address: CERN, CH-1211 Geneva 23, Switzerland. 

and 

1 Dipartimento di Fisica dell'Universita di Bari, Bari, Italy. 

S INFN Sezione di Bari, Bari, Italy. 

1 Dipartimento di Fisica dell'Universita della Calabria, Cosenza, Italy. 

U INFN Gruppo collegato di Cosenza, Cosenza, Italy. 

v Department of Nuclear and Particle Physics, Uppsala Univeristy, Uppsala, 

Sweden. 

w A. Soltan Institute for Nuclear Studies, Warsaw, Poland. 



Abstract 

We report the first observation of the rare n — > e + e~e + e~ decay based on 1.7 fb _1 
collected by the KLOE experiment at the DA<3?NE ^-factory. The selection of the 
e + e~e + e~ final state is fully inclusive of radiation. We have identified 362 ± 29 
events resulting in a branching ratio of (2.4 ± 0.2 sta t+bckg ± 0A syst ) x 10~ 5 . 
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1 Introduction 



The T] — > e + e~e + e~ decay proceeds through two virtual photons intermediate 
state with internal photon conversion to e + e~ pairs. Conversion decays offer 
the possibility to precisely measure the virtual photon 4-momentum via the 
invariant mass of the e + e~ pair. The lack of hadrons among the decay prod- 
ucts makes the matrix element directly sensitive to the 7] meson transition 
form factor [TJ. The knowledge of the rj meson coupling to virtual photons 
is important for the calculation of the anomalous magnetic moment of the 
muon, being pseudoscalar exchange the major contribution to the hadronic 
light-by- light scattering. 

The first theoretical evaluation, T(r] — > e + e~e + e~)/r(r] — > 77) = 6.6 x 10~ 5 , 
dates from 1967 [2j. The width ratio translates into a branching ratio (BR) 
BRij] — > e + e~e + e~) = 2.59 x 10~ 5 when the world average of the BR(i] — > 77) 
measurements [5] is taken as normalization factor. Other predictions exist in 
literature [41I51I61IT] . with differences at the level of 10%. 

Double lepton-antilepton rj decays have been searched by the CMD-2 and 
the WASA experiments, obtaining the upper limits at 90% C.L., BR(r] — > 
e + e~e + e-) < 6.9 x 10~ 5 [8] and BR(rj -> e + e~e + e-) < 9.7 x 10" 5 [9j, respec- 
tively. 



2 The KLOE detector 



The KLOE experiment operates at DA$NE , the Frascati 0-factory. DA$NE 
is an e + e~ collider running at a center of mass energy of ~ 1020 MeV, the 
mass of the 4> meson. Equal energy positron and electron beams collide at an 
angle of 7T-25 mrad, producing nearly at rest </> mesons. 

The detector consists of a large cylindrical Drift Chamber (DC), surrounded by 
a lead-scintillating fiber electromagnetic calorimeter. A superconducting coil 
around the EMC provides a 0.52 T field. The beam pipe at the interaction 
region is spherical in shape with 10 cm radius, it is made of a Beryllium- 
Aluminum alloy of 0.5 mm thickness. Low beta quadrupoles are located at 
about ±50 cm distance from the interaction region. The drift chamber [TO"] . 
4 m in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense wires and 
37,746 aluminum field wires. The chamber shell is made of carbon fiber-epoxy 
composite with an internal wall of 1.1 mm thickness, the gas used is a 90% 
helium, 10% isobutane mixture. The spatial resolutions are a xy ~ 150 /zm and 
a z ~ 2 mm. The momentum resolution is cr{p±)/p± ~ 0.4%. Vertices are recon- 
structed with a spatial resolution of ~ 3 mm. The calorimeter [TTJ is divided 
into a barrel and two endcaps, for a total of 88 modules, and covers 98% of the 
solid angle. The modules are read out at both ends by photomultipliers, both 
in amplitude and time. The readout granularity is ~ (4.4 x 4.4) cm 2 , for a to- 
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tal of 2440 cells arranged in five layers. The energy deposits are obtained from 
the signal amplitude while the arrival times and the particles positions are ob- 
tained from the time differences. Cells close in time and space are grouped into 
calorimeter clusters. The cluster energy E is the sum of the cell energies. The 
cluster time T and position R are energy-weighted averages. Energy and time 
resolutions are a E /E = 5.7%/ \j E (GeV) and o t = 57 ps/^E (GeV) © 100 ps, 
respectively. The trigger [12] uses both calorimeter and chamber information. 
In this analysis the events are selected by the calorimeter trigger, requiring 
two energy deposits with E > 50 MeV for the barrel and E > 150 MeV for 
the endcaps. A cosmic veto rejects events with at least two energy deposits 
above 30 MeV in the outermost calorimeter layer. Data are then analyzed by 
an event classification filter [T3], which selects and streams various categories 
of events in different output files. 



3 Event selection 

The analysis has been performed using 1,733 pb _1 from the 2004-2005 data 
set at y/s ~ 1.02 GeV. 242 pb _1 of data taken off-peak at y/s = 1.0 GeV 
were used to study the e + e~ continuum. Monte Carlo (MC) events are used 
to simulate the signal and the background. The signal is generated accord- 
ing to the matrix element in [5], assuming BR = 2.7 x 10~ 5 , in a sample of 
167,531 pb _1 . Other MC samples are: 3,447 pb _1 simulating the main (p de- 
cays (<ft — > KK and <fi — > pir) and 17,517 pb -1 simulating others more rare 
<j) decays. All MC productions account for run by run variations of the main 
data-taking parameters such as background conditions, detector response and 
beam configuration. Data-MC corrections for calorimeter cluster energies and 
tracking efficiency, evaluated with radiative Bhabha events and <fi — > pit sam- 
ples respectively, have been applied. Effects of Final State Radiation (FSR) 
have been taken into account using the PHOTOS MC package [T4]fT5] . This 
package simulates the emission of FSR photons by any of the decay products 
taking also into account the interference between different diagrams. PHO- 
TOS is used in the Monte Carlo at the event generation level, so that our 
simulation fully accounts for radiative effects. 

At KLOE, 7] mesons are produced together with a monochromatic recoil pho- 
ton (E^ = 363 MeV) through the radiative decay <p ~^ VI- I n the considered 
data sample about 72 x 10 6 rfs are produced. As first step of the analysis, a pre- 
selection is performed requiring at least four (two positive and two negative) 
tracks extrapolated inside a fiducial volume defined by a cylinder centered in 
the interaction point and having radius R = 4 cm and length Az = 20 cm. 
For each charge, the two tracks with the highest momenta are selected. One 
and only one neutral cluster, having energy E c i > 250 MeV and polar angle in 
the range (23° — 157°), is required. A cluster is defined neutral if it does not 
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have any associated track and has a time of flight compatible with the photon 
hypothesis. To improve the energy and momentum resolution, a kinematic fit 
is performed imposing the four-momentum conservation and the photon time 
of flight. A very loose cut on the \ 2 °f the kinematic fit (x 2 < 4000) is applied 
in order to discard poorly reconstructed events. 



4 Background rejection 

Two sources of background are present: 

(1) <p background: 

this is mainly due to cf) — > 7r + 7r~7r° events (with 7r° Dalitz decay) and 
to — > 777 events either with rj — > 7r + 7r~7r° (with tt° Dalitz decay) or 
f] — > 7r + 7r~e + e~ or with rj — > e + e~7 (with photon conversion on the Beam 
Pipe, BP, or the DC inner Wall, DCW). This last background has the 
same signature of the signal. Background from — > KK is also present 
at the preselection level. 

(2) e + e~ continuum background: 

this is mainly due to e + e~ — > e + e^( , ~f) events with photon conversions, 
split tracks or interactions in the DA$NE low beta quadrupoles. This 
background has been studied using off-peak data taken at = 1 GeV, 
where <fi decays are negligible. 

A first background rejection is performed cutting on the sum of the absolute 
value of the momenta of the four selected tracks requiring (600 < J2i \Pi\ < 
700) MeV. 

To remove e + e~ continuum background from interactions in the low beta 
quadrupoles, the quantities (cos9f) and (cos 9b) have been defined as the av- 
erage polar angle of forward and backward selected particles. Events having 
(cos 6f) > 0.85 and (cos 6b) < —0.85 are rejected. This cut has no effect on 
signal selection efficiency. 

To reject events due to photon conversion, each track is extrapolated backward 
to the intersection with the BP and with the DCW. For each track pair, the 
invariant mass (M ee ) and the relative distance (D ee ) are computed. A clear 
signal of photon conversion is visible in the D ee -M ee 2D plot for BP and DCW 
(figured]). Events having at least one combination satisfying M ee (BP) < 10 
MeV and D ee (BP) < 2 cm or M ee (DCW) < 30 MeV and D ee (DCW) < 2 cm 
are rejected. 

The last rejection is based on the Particle IDentification (PID) of charged par- 
ticles. For each track associated to a calorimeter cluster, the quantity At = 
Urack—tduster in both electron and pion hypothesis is evaluated; t trac k is defined 
as the length of the track divided by /3(m)c. Track with At e / At n < 1(> 1) 
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Table 1 

Number of events in data, MC signal efficiency, background rejection factor at dif- 
ferent steps of the analysis. 



Cut 


N(data) 


F~.* n ( SI P* 1 
°ana \ OL SJ 


R + - 


R<j)-+KK/ fm 


-^others 


Pvpcplpptinri 

J. I CoClCU IHJ11 


4^1 Q94 


9R C SM ^ 


1 Sfi^9"i v in 2 


^ ni o\ x 1 n 3 


1 d^fR) Y ID 3 


x 2 


36282 


0.217(1) 


2.01(3) x 10 2 


1.13(1) x 10 5 


3.44(5) x 10 4 




16811 


0.216(1) 


2.68(5) x 10 2 


2.21(3) x 10 5 


6.9(1) x 10 4 


cos#t,, cos 6 f 


15003 


0.216(1) 


2.68(5) x 10 2 


2.21(3) x 10 5 


6.9(1) x 10 4 


7 conversion 


12198 


0.209(1) 


1.11(4) x 10 3 


2.53(3) x 10 5 


1.13(3) x 10 5 


PID 


4239 


0.205(1) 


1.12(4) x 10 3 


1.02(8) x 10 7 


5.1(3) x 10 5 



are identified as electron (pion). Events having more than two pions or no 
electrons are discarded. 

The effects of background rejection cuts on the various data components are 
visible in figure EJ where the four electrons invariant mass, M e + e - e + e - , is 
shown at different steps of the analysis. In table HJ number of events in data, 
N(data), MC signal efficiencies, £ ana (sig), and background rejection factor R, 
defined as the ratio of analysis efficiency between signal and background, are 
also reported. The R value has been evaluated for three different cathegories: 
->■ 777 with V -> e + e"7 (i^ e+e - 7 ), -> KK and -> p?r (R^kR/^) 
and all other decays products (i2 D thers)- After all cuts, background from 
kaons and — > 7r + 7r~7T° events is negligible. The same holds for all other 
decays but 77 — > e + e~7 which, as will be shown in the next section, results 
in ~ 15% contamination level. Systematics on the Monte Carlo description of 
photon conversion have been studied using events with similar characteristics. 
A clean control sample is provided by the — > r]e + e~, 77 — > 7r + 7r~7r° decay 
chain, where simple analysis cuts provide a good data-MC agreement, with 
negligible background contamination. As for the 77 — > e + e~e + e~ channel, be- 
fore dedicated analysis cuts the control sample is significantly contaminated 
by background from photon conversion (0 — > 777 with photon converting on 
beam pipe or drift chamber walls). This background is completely removed 
rejecting events with D ee (DCW) < 10 cm and M ee (DCW) < 80 MeV. For 
the 77 — > e + e~e + e~ channel this cut has not been applied because, having 
two electrons and two positrons in the final state, the search for a conversion 
has to be performed over all the four e + e~ combinations, thus spreading the 
signal contribution in the M ee (DCW)-D ee (DCW) plane and lowering signif- 
icantly the analysis efficiency. Removing the cuts on M ee -D ee planes in the 
control sample, a clear background contamination from photon conversion 
is visible. Data-MC comparison shows that, increasing in the simulation the 
probability of conversion by 10%, an excellent agreement is found. A 10% 
systematic error is then assigned to photon conversion and added to the un- 
certainties coming from MC statistics and BR(t7 — > e + e~ r f) measurement [3J: 
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Fig. 1. D ee vs M ee evaluated at the drift chamber wall for MC — > rjj 
background (top panel) and MC signal (bottom panel). Events in the box 
M ee (DCW) < 30 MeV n D ee (DCW) < 2 cm are rejected. 



N( V -> e+e'y) = 80 ± 3 MC ± 8 BR ± 8 



syst- 



5 Evaluation of the BR(t7 — > e+e e+e (7)) 



As discussed in the previous section, the only significant background con- 
tamination surviving all the analysis cuts is due to r\ — > e + e~7 events with 
photon conversion, that have a signature similar to the signal. The overall 
estimated background from decays has been subtracted bin-by-bin to the 
M e + e - e + e - spectrum obtained in data (figure [3] top), taking into account also 
systematic errors. The event counting is done fitting the resulting spectrum 
with the two residual contributions: signal and e + e~ continuum background 
events. The M e + e - e + e - shape for the signal is obtained by fitting MC events 
with two Gaussian functions plus a third order polynomial function. The fit 
range is 500 < M e + e - e + e - < 600 MeV. The M e + e - e + e - distribution for e + e~ 
continuum events has been studied on the data taken at y/s = 1 GeV, where 
contributions from <ft decays are suppressed. Even though the small statis- 
tics of the sample does not allow to precisely extract the shape, a first order 
polynomial well reproduces the data in the signal region. The free parameters 
are an overall scale factor for signal and the two parameters describing the 
e + e~ continuum background. Fit results are shown in figure [3] bottom. The 
resulting % 2 /ndf is 43.9/34, corresponding to P(% 2 ) = 0.12. The number of 
signal events is N(rj — > e + e~e + e~) = 362 ± 29. The branching ratio has been 
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M [MeV] M [MeV] 

eeee L J eeee L J 

Fig. 2. M e + e - e + e - distribution after different analysis cuts: white: after the Wi\ 
and the {cos 9) cuts; gray: after the cut on photon conversion; black: after the PID 
requirement. Top left: data; top right: off-peak; bottom left: <fi background Monte 
Carlo; bottom right: signal Monte Carlo. 

evaluated according to the formula: 

B R(r, -> e + e-e+e-( 7 )) = jVgtpt^W 1 (1) 

where N v ^. e + e - e + e -^ is the number of signal events and e v ^. e + e - e + e -(^ is the 
efficiency taken from MC. The number of — > 777 events, N V1 , has been ob- 
tained using the formula N V1 = C-a^ rni where £ is the integrated luminosity 
and the cross section a^ vl has been evaluated taking into account the me- 
son line shape on a run by run basis [TH] . Inserting all the numbers quoted in 
table [21 the value: 

BR(r, -> e + e-e + e-( 7 )) = (2.44 ± 0.19 stat+bckg ) x 1(T 5 (2) 

is obtained, where the error accounts for the uncertainty of the fit result. 
The systematic uncertainties due to analysis cuts have been evaluated by ap- 
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Fig. 3. Top panel: M e + e - e + e - data distribution at the end of the analysis chain; the 
expected (ft background MC shape is shown in dark gray. Bottom panel: M e 
fit to data after (ft background subtraction. 



Table 2 

Summary of the numbers used in the master formula ([1]) for the branching ratio 
evaluation. 



BR inputs 



Values 



Number of events 
Efficiency e v ^ e + e - e + e -( 7 ) 
Luminosity 



362 ± 29 
0.205 ± 0.001 
^1733 ± 10) pb" 1 



e^e 



rpi cross section (41.7 ± 0.6) nb 



plying separately a variation of ±1 a on all variables and re-evaluating the 
branching ratio. The o values have been obtained using MC signal events. For 
the x 2 variable the cut has been moved by ±500, while the particle identifi- 
cation cut has been changed by ±10%. The systematic error on the fit to the 
M e + e - e + e - distribution has been evaluated considering: 

• the binning of the M e + e - e + e - histogram, changed from 3 MeV, used as 
default, to 2 and 4 MeV; 

• the M e + e - e + e - range, enlarged and reduced by 10 MeV on both sides; 

• the slope of the e + e~ continuum background has been fixed to the value 
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Table 3 

Summary table of systematic uncertainties. 



Source of uncertainty 


Relative error 


x 2 


-0.51% / + 2.62% 


(cos 6b) and (cosOf) 


-0.04% / + 0.47% 


Ei \Pi\ 


+0.11% 


7 conversion 


-0.74% / + 2.43% 


PID 


+1.84% 


Fit range 


-0.38% / + 1.13% 


Binning on M e + e - e + e - 


-3.21% / + 0.19% 


Background slope 


+0.38% 


Normalization 


±1.64% 


Total 


-3.73% / + 4.53% 



obtained from off-peak data fit. 

The relative variation of the BR for each source of systematic uncertainty is 
reported in table [3j The uncertainty on N V1 has been added to the systematics 
in the normalization term. The total error is taken as the quadratic sum of all 
contributions. 



6 Conclusions 



Using a sample of 1.7 fb collected in the (j> meson mass region, the first 
observation of the r\ — > e + e~e + e~ (7) decay has been obtained on the basis of 
362 ± 29 events. The corresponding branching ratio is: 

BR( V -+ e + e-e + e-( 7 )) = (2.4 ± 0.2 stat+bckg ± 0.1 syst ) x 1(T 5 . (3) 

Radiative events slightly modify momentum distribution of the charged par- 
ticles and have been carefully considered in the efficiency evaluation. As a 
result, the measured branching ratio is fully radiation inclusive. 
Our measurement is in agreement with theoretical predictions, which are in 
the range (2.41 - 2.67) x 10" 5 [214115116117] . 
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